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We study mode-locking in disordered media as a boundary-value problem. Focusing on the sim- 
plest class of mode-locking models which consists of a single driven overdamped degree-of-freedom, 
we develop an analytical method to obtain the shape of the Arnol'd tongues in the regime of low 
ac-driving amplitude or high ac-driving frequency. The method is exact for a scalloped pinning po- 
tential and easily adapted to other pinning potentials. It is complementary to the analysis based on 
the well-known Shapiro's argument that holds in the perturbative regime of large driving amplitudes 
or low driving frequency, where the effect of pinning is weak. 

PACS numbers: 



I. INTRODUCTION 

The phenomenon of mode-locking is a general feature 
of nonlinear dynamical systems. It consists of the res- 
onant response to an external periodic force that oc- 
curs when a characteristic frequency of the driven sys- 
tem matches or locks onto the driving frequency. In the 
mode-locked region the system traces periodic orbits in 
phase space. Outside the region of mode locking, the 
system may follow quasiperiodic orbits or march towards 
the onset of chaos. 

Examples of systems that exhibit mode-locking 
abound in nature. In 1665, Huygens discovered the spon- 
taneous synchronization of swinging pendulum clocks in 
close proximity to one another. On the celestial scale, the 
moon's period of rotation locks onto its period of revo- 
lution about Earth in a l:l-ratio, so it always presents 
the same face to Earth. The rotation of Mercury is also 
locked onto the Sun such that there are three rotations for 
every two orbits, constituting a 3:2-mode locking. Mode- 
locking to external periodic stumuli is a common feature 
in biology. Examples include the cell cycle in budding 
yeast [1], the swimming and heartbeat networks of medic- 
inal leeches [5], and the rhythmic behavior produced by 
neuronal networks [5]. Mode-locking is also relevant in 
such diverse phenomena as vortex shedding [4 , singing 
sand dunes [5], and multimode lasers [B]. In condensed- 
matter physics, Josephson junction arrays [7], driven su- 
perconducting vortices [H H] , and charged-density waves 
(CDW) m [H Uni [HI [15] provide convenient settings 
for studying mode-locking. 

Mathematicians have long used the language of maps 
to describe the trajectories followed by dynamical sys- 
tems in phase space and thus infer their physical prop- 
erties. In particular, circle maps ^fWl [T5] . along with 
the mathematical tools of bifurcation theory and return 
map, lay the foundation for much of the theoretical mod- 
eling and analysis of mode-locking phenomena. Concep- 
tually, the simplest evolution equation that yields mode- 



locking consists of a single overdamped degree-of-freedom 
(DOF) in a periodic pinning potential Vp^cj)) of stength 
h, driven by an external periodic force of frequency uj. 
The dynamics is described by the equation 

^ = Fo + Fi cos cot + Fp{cp) , (1) 

where Fp{(j)) = —h"^^^^ is the pinning force. Mode- 
locking occurs when the time-averaged velocity locks on 
to a rational multiple of the external driving frequency, 

/d(j)\ f 1 d(t> 

X'^^ /cycle" Avc^ 27r/cj 

= Fo+{Fp)^^^,^^u;,= ^u; (2) 

for a region of nonzero area in the {Fq, _F'i)-space. The re- 
gions of the {Fq, Fi) parameter space where mode-locking 
occurs are known as Arnol'd tongues [T7]- The time- 
averaged velocity exhibits a devil's-staircase (DS) behav- 
ior: there exists a mode-locked plateau corresponding to 
each rational p/q. While there have been much work 
on the fractal dimensionality of the set of gaps between 
the mode-locked steps in the staircase [TH] [THl [10] , the 
calculation of the width of the mode-locked steps has 
mainly relied on an argument put forward by Shapiro for 
a single-particle model of CD Ws |2T] j^l] . In this case the 
DOF (j) is the phase of the CDW, while Fq and Fi are 
the amplitudes of the dc and the ac driving electric fields, 
respectively. The time-averaged velocity corresponds to 
the drift velocity cu^ of the CDW, which in turn deter- 
mines the CDW current. 

Most available results on mode locking have been ob- 
tained numerically. One of the rare analytical results is 
for a single DOF in a cosine pinning potential, Vp{(f>) — 
h cos (j). In this case the Arnol'd tongues are symmetric 
with respect to the mirror axis centered at the apex of 
the tongue and pinch to zero width for all parameters in 
its phase space, as shown in Fig. [T] On the other hand. 



2 



numerics have shown that in extended systems consisting 
of many coupled DOFs the Arnol'd tongues are generally 
asymmetric and never pinch to zero width. In this work, 
we reexamine single DOF model for various pinning po- 
tentials and show that the cosine pinning is special. For 
a generic pining potential the tongues are asymmetric 
and do not pinch to zero as soon as the pinning poten- 
tial contains any higher harmonics. This more complex 
shape of the Arnol'd tongues may be the result of the 
collective behavior of many coupled degrees of freedom, 
but it arises even for a single particle, provided the pin- 
ning potential differs from a simple cosine. Of course in 
an extended system collective effects renormalize the pin- 
ning potential, so that even a bare cosine pinning yields 
asymmetric Arnol'd tongues. We show below that for 
generic pinning potential the Shapiro's method always 
fails at small Fi, where the asymmetry of the Arnol'd 
tongue is most apparent. For the specific case of scal- 
loped potential, we develop an exact analytical mehtod 
for the calculation of the mode-locking behavior of one 
DOF for Fi small enough that the single DOF does not 
hop from one scallop to the next. This is precisely the 
regime where the Shapiro's argument always fails. Our 
analytical solution fro small Fi also provides dynamical 
constraints that are both necessary and sufficient for the 
full determination of the mode-locked step widths. 

This work serves as a starting point in our effort to con- 
struct a mean-field theory for the general phenomenon of 
mode-locking in extended media that is amenable to an- 
alytical analysis in useful limits. To do so, we must first 
test analytical approximations at the single-particle level 
and identify a simple, yet generic, pinning potential as 
our prototypical starting model. To achieve these two 
first-step goals, we study in detail the case of l:l-mode 
locking and consider two specific periodic pinning po- 
tential G'[0], a scalloped parabolic and an impure cosine 
potential. In section II we review Shapiro's method for 
calculating the width of the mode-locked steps. We then 
consider the one-particle mode-locking dynamics in the 
scalloped parabolic pinning potential in section III. This 
form of potential has the advantage of permitting exact 
analytical solution in the no-hopping regime. In section 
IV, we repeat our analysis for the impure cosine pinning 
potential. Section V concludes our paper. 



II. SHAPIRO'S METHOD IN CALCULATING 
MODE-LOCKED STEP WIDTHS 

Adapting an argument proposed by Shapiro |21j . 
Thorne et. al. computed the width of the harmonic and 
subharmonic mode-locked steps for the one-DOF model 
of CDW given in Eq. ^ [25] . While their calculation ex- 
plains the occurrence of mode- locking in terms of the low- 
ering of the pinning energy, their results do not agree with 
numerical results for the shape of the Arnol'd tongues in 
the regime of low ac-driving amplitude or high ac-driving 
frequency. This regime corresponds to the case where the 



particle dynamics is confined to one period of the pin- 
ning potential. To show how Shapiro's method fails in 
this regime, we first review the calculation of Thorne et 
al. [22^. 

It is instructive to first consider the single DOF model 
of Eq. (jl} for a cosine pinning potential. We note that 
in the absence of pinning the equation of motion has the 
exact solution (/){t) — FqI + — sinwi + (j)o, with (/)q a 
constant, which gives the trivial result LUd = Fq. In the 
presence of pinning we let 



(f){t) = UJdt - 



Fi 

sin LUt + 00 (^) : 



(3) 



with ujd and 4'o{t) to be determined. We note that 
Eq. (jlj contains two characteristic frequency scales: the 
frequency to of the external drive and the frequency h 
of temporal variations of the phase due to the pinning 
potential. In the limit of large drive, the particle moves 
rapidly over the pinning potential and the temporal vari- 
ations due to pinning are small. Following Ref. [55] we 
then look for a solution where 0o ~ constant, indepen- 
dent of time. We expect this approximation will apply for 
large Fi and weak pinning. It turns out to be essentially 
exact for a simple cosine pinning potential, but it fails for 
small values of Fi for arbitrary periodic pinning. Substi- 
tuting Eq. ([3| into Eq. (jlj, utilizing the Besscl function 
summation formula 



exp 



sin u!t 



p— — oo 



exp {—ipujt) , 



(4) 



and averaging over a cycle, we obtain 



u:d^F,+ {F,)^^^,^ , (5) 

OO 



p=i 



where the sum is over all p such that LOd — pto. Similarly 
the mean pinning energy is given by 



cycle 



hcos ( 



oo 

E 

p=l 

(-d=p-) 



i-irj,iF,/Lo) . (7) 



A mode-locked state occurs when ujd in Eq. ([5| can be 
kept constant (and equal to puj) for a range of values of 
_Fo and Fi by adjusting the phase 0o- For a cosine pin- 
ning potential the range of c/jq that renders tOd — constant 
is identical to the range of 0o where {Vp}^^^^^ < 0. In this 
case the region of parameters where the system is mode- 
locked coincides with the region where the energy of the 
mode-locked state is lower than that of the unlocked 
state. This is not, however, the case for arbitrary pinning 
potential, as we will see below. In general the condition 
that energy of the mode-locked state be lower than that 
of the unlocked state is necessary for their stability of 
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FIG. 1: 1:1 Arnold tongue for pure cosine pinning potential. 



the mode-locked state, but not sufficient (see Fig. [5] be- 
low). For pure cosine pinning case only harmonic steps 
are obtained. The p/1 step occurs for a range (AFo)p/i 
of values of Fq given by (Ai^o)p/i = 2/i| Jp(Fi/w)|. The 
Arnol'd tongues for this case are shown in Fig. [l] 

Thorne and collaborators adapted this argument to 
an arbitrary periodic pinning potential that can be ex- 
pressed as a Fourier sum: 



9=1 



(8) 



Proceeding as in the single cosine case, using again the 
Bessel summation formula and averaging over a cycle, 
we find that uj^ is determined by Eq. (|5|, where the 
time-averaged pinning force is zero unless LUd = P^/q, 
for which 



cycle 



oo 
9=1 P 



(9) 



Thus, for a given Fi, the constant phase 0o can be ad- 
justed to compensate changes in Fq such that the time- 
average phase velocity ujd stays locked to rational val- 
ues of the external driving frequency. Furthermore, Eq. 
^ suggests that the width of the Shapiro steps oscil- 
lates with Fi, which has indeed been observed in exper- 
iments [HI [TU] . These authors further assumed that the 
range of parameters where the system is mode-locked co- 
incides with those where the mean pinning energy in the 
mode locked state is lower than its value ao/2 in the un- 
locked state, i.e.. 



oo oo 



S{Vp)cycie = X! X! aqjp COS {q4>o) <0 

q—1 P— — 00 



(10) 



The mean pinning force given by Eq. ([9| cancels the 
DC drive Fq in a range —(j)p{Fo, Fi) < (f>o < (t)p{Fo,Fi). 



FIG. 2: The scalloped parabolic pinning potential, defined by 
Eq. pi. 



The condition that 0p(-F'o, Fi) e [0, 2tt] can then be used 
to determine the boundaries of the mode-locked regions 
in the (Fo,Fi) plane. Thorne et al. use the condition 
6{Vp)cycie < to determine the Arnold tongues. The lat- 
ter holds for a range -4im{Fo,Fi) < (f>Q < (j),n{Fo, Fi). 
For pure cosine pinning (f>„^{FQ,Fi) — (j)p(FQ, Fi). For 
arbitrary pinning potentials 4>m{FQ, Fi) > 4>p{Fq,Fi), re- 
sulting in an overestimate of the mode-locked regions. 
Regardless of the condition used for mode-locking, a con- 
sequence of Eqs. (|9| and ( 10 1 is that each Arnol'd tongue 
is symmetric about a central axis for any periodic pinning 
potential. In the case of l:l-mode locking, the tongue is 
symmetric about the axis Fg = 1, since Jp(0) = for all 
p ^ and Jo(0) = 1 assuming appropriate normalization. 
For the l:l-mode locking case, we will show explicitly, 
via exact numerical calculations, that this symmetry of 
Arnol'd tongue is violated in the regime where the influ- 
ence of ac-drive is large, i. e. when the drive amplitude 
Fi is small or the drive frequency uj is small (c./. Fig. |4|, 
for two specific cases of the periodic pinning potential 
corresponding to the scalloped parabolic and the impure 
cosine pinning. 



III. SCALLOPED PARABOLIC PINNING 
POTENTIAL 

We first consider the case of scalloped parabolic pin- 
ning, whose potential can be expressed either in terms of 
the floor function, [x\ , defined as the largest integer less 
than or equal to its argument x, or as a Fourier series: 



Kicallopod(</') 



L0J 



24 ^ 27r2fl2 



cos 27rg0 



(11) 
(12) 



A plot of the scalloped parabolic potential is shown in 
Fig. |2] In this case the double sum in the expression for 
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the time-averaged pinning force, Eq. ([9]), reduces to a 
single infinite sum over terms for which p = q. Using Eq. 



(12 1, the cycle-averaged pinning force is given by: 



{Fp{<f>o)) cycle = — Jq(gFi) sin27r#o 



(13) 



q=Q 



For the scalloped pinning potential of Fig. |2]the pinning 
force is piece-wise linear and the dynamics of the driven 
particle can be solved exactly within each period of the 
pinning potential. This solution is expected to be exact 
for small values of Fi, provided the particle does not 
hop from one scallop to the next i the time scale of the 
external drive. Within one period Eq. ([T|) is simply 



dt 



Fn 



Fi cos 27ri 



(14) 



For simplicity we discuss in detail only mode-locked steps 
with p = 1. In this case Eq. ( 14 1 must be solved with the 
boundary conditions 



<t){tj + nT) ^ 
(f>{tj + {n+l)T) = 1 



(15) 
(16) 



where T — 2Trq/uj. The jump time tj plays the role of 
the constant phase (j)Q in the previous section. Equa- 
tion ( 15 ) determines the constant of integration for our 

( 16 1 yields the relationship 
This 



first-order ODE, while Eq. 



between Fq and -Fi corresponding to mode-locking 
can be rewritten as 



1 



Fo 



Fi 



2 ^cos 2ntj 



27rsin 2Trtj) 



The condition that the values of tj in Eq. (17 1 must be 



in [0, T] determines the boundaries Fi{Fq) of the Arnol'd 
tongues, given by 



Fi > ±vTT(2^ 



(18) 



These are the straight lines shown in Fig. [3] along with 
the complete Arnol'd tongues for the 1:1 step obtained by 
exact numerics. Clearly the analytical solution yields the 
exact value of Fq for the onset of mode-locking at Fi = 
and also does an excellent job of fitting the exact Arnol'd 
tongues for small Fi, where the driven particle remains 
within a single scallop. In contrast, Shapiro's argument 
fails most severely precisely in this region of small Fi, 
as apparent from Fig. [4] where the Arnol'd tongues ob- 
tained for the 1:1 mode-locked step by the Shapiro argu- 
ment (triangles) are compared to the exact solution (di- 
amonds). As mentioned in section II, Shapiro's method 
predicts symmetrical oscillation about the axis Fq — 1. 
However, exact numerical solution reveals that the 1:1- 
mode locked step, in fact, originates from Fq « 1.1. This 
difference arises because the Shapiro method is essen- 
tially a perturbation theory about the high velocity state 




FIG. 3: 1:1 Arnold tongue for scalloped parabolic pinning 
potential: Thome's method versus "exact" method for small 



(17) 




FIG. 4: Plot of Arnol'd tongue for 1:1 mode-locking step with 
the scalloped parabolic pinning potential. Black triangles rep- 
resent data points obtained from the zeros of the pinning en- 
ergy using Thome's method. Numerically, the first 100 terms 
of the sum were included in our summation over the Bessel 
functions. The black diamonds represent data from numerical 
results. 



and becomes exact at large drives and very weak disorder. 
On the other hand for fixed driving frequency u and small 
ac-drive amplitude Fi, or conversely fixed Fi and high 
frequency the driving force has only a weak effect and 
pinning dominates. In this region the time-averaged ve- 
locity is well approximated by the instantaneous velocity. 
This fact underlies our approximation scheme in which 
we solve exactly for the instantaneous dynamical solu- 
tion for the CDW phase (j>{t), assuming that the particle 
does not hop over any one scallop in the periodic pinning 
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FIG. 5: Plot of the impure cosine pinning potential consisting 
of a sum of two harmonics. Gray triangles represent data 
points obtained from the zeros of the pinning energy using 
Thome's method. The black diamonds represent data from 
numerics. 



FIG. 6: Shape of the Arnol'd tongue for the 1:1 mode-locking 
step for an impure cosine pinning potential. The diamonds 
are exact numerical results and the V-shaped straight lines 
are the analytical approximation in the no-hopping regime, 
corresponding to the limit of small drive/high frequency. 



potential. Neglecting the effect of scallop-hopping simpli- 
fies the dynamics considerably, as the elimination of the 
associated nonlinearity permits analytically tractable so- 
lutions. This method is complementary to the one devel- 
oped by Shapiro. Our method can be readily generalized 
to other harmonic and subharmonic mode-locking steps 
that involve no hopping between scallops and thus satisfy 
the constraint p = 1. 



IV. IMPURE COSINE PINNING POTENTIAL 

For a cosine pinning potential the width of mode- 
locked steps pinches to zero periodically at finite values 
of the ac-driving amplitude, Fi. This is well explained 
by Shapiro's argument [22] based on its expression of the 
time-averaged pinning potential: since the time-averaged 
pinning force has only one term for any particular values 
of p and q and it is a consequence of the fact that pinning 
potential contains only a single harmonic in its Fourier 
series. To see this point explicitly, we compare the two 
cases of a pure cosine pinning potential function and of 
an impure cosine, consisting of a sum of two different 
harmonics: 

Vcosino('/>) = -!-cos27r0, (19) 

ZTT 

cos27r0+ cos47r(/) . (20) 



(</>) 



impure \ T/ r, -r ' A 

Again we are focusing on l:l-mode locked steps. For the 
impure cosine pinning, the time-averaged pinning force 
is 



{FpiM) Ji(^^i)sin27r0o 



-O.IJ2 (2Fi)sin47r0o 

(21) 



It is apparent from Fig.[5]that the addition of a small har- 
monic is sufficient to give a finite width to the tongues 
for all nonzero values of Fi. Again Shapiro's method 
clearly fails at small Fi. To analyze the dynamics in 
this region, we assume that there is no particle-hopping 
and use the analytical method discussed in the previ- 
ous section. To do this we fit the impure cosine pinning 
potential to a parabola over one cycle. Using Mathe- 
matica, we find that the fitted pinning potential yields 
Kosino(<?^) = 0.186-1. 170-1-1. 19(/)2, or equivalently, a scal- 
loped parabolic pinning potential with a pinning strength 
of /i « 2.39 over the same cycle. The resulting com- 
parison plot of the numerically obtained Arnol'd tongue 
for l:l-mode locking and the "V" -shaped asymptotes ob- 
tained analytically in the no-hopping approximation is 
shown in Fig. (|6]). Again, there is remarkable agreement 
between our analytical approximation and the exact nu- 
merical solution in the small-drive regime. 

To further clarify the role of the no-hopping approxi- 
mation used in our analytical calculation, we note that in 
obtaining Eq. (18 1, we implicitly assume that \4>{t)\ < 1, 
for t g [0, 27r]. Self-consistency of our solution is then 
imposed by looking for the set of points in the region 
defined by Eq. (18 1 that explicitly satisfy \(j>{t)\ < 1 for 
t £ [0,2n]. As shown in Fig. [t] the region of {Fo,Fi), 
bounded by black circles, which explicitly satisfies the 
constraint \4>{t)\ < 1 is an excellent approximation to 
the lopsided pear-shaped region of the Arnol'd-tongue 
for Fi < 4. Beyond this region, the particle presumably 
hops over more than one scallop per drive period, and 
our analytical approximation would no longer apply. 
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Fi/h 1:1 




0.6 0.8 1 1.2 1.4 

Fo/h 

FIG. 7: The shape of the Arnol'd tongue for the 1:1 mode- 
locking step for a scalloped paraboHc pinning potential. The 
solid line are exact numerical results and the (red onHne) V- 
shaped straight lines correspond to the analytical approxima- 
tion in the no hopping regime. The region bounded by the 
black circles are points of (-Fo, ^^i) which are exphcitly checked 
to satisfy the constraint of \(f>{t)\ < 1. 



V. CONCLUSION 

In summary, we have developed a simple analytical 
method to obtain the shape of the Arnol'd tongues in the 
regime of small ac-driving amplitude Fi or high-driving 



frequency, where the driven particle does not hop be- 
tween different periods of the driving potential within one 
period of the drive. This method is complementary to the 
perturbative one based on Shapiro's argument f5TJ [52] 
that applies in the large Fi or low frequency regime. The 
method is exact for a scalloped pinning potential and is 
easily adapted to other pinning potentials by a simple 
fit. Our method is easily adapted to the analysis oip/q 
mode-locking steps for arbitrary q and p ~ I- 

As mentioned in the Introduction, the motivation of 
this work was to develop simple methods for the analy- 
sis of mode-locking steps that will serve as the starting 
point for the development of a mean-field theory of mode- 
locking in systems composed of many interacting many 
degrees-of-freedom. We hope to do this by combining 
the no-hopping approximation for the analysis of the low 
drive regime with the Shapiro argument for the study 
of the high-drive region. Finally, while the phenomenon 
of mode-locking is pervasive in dynamical systems, each 
subfield has developed its own set of theoretical tools, 
often inspired and required by real-world applications. 
Despite past attempts a consensus is still lacking 
in regard to a unifying formalism in the description of 
mode-locking. It is our further hope that our results 
would illuminate aspects of mode-locking, which are still 
not well understood and which would merit further study. 
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